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High resolution AMI Large Array imaging of spinning dust sources: 
spatially correlated 8 fim emission and evidence of a stellar wind in 
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ABSTRACT 

We present 25" resolution radio images of five Lynds Dark Nebulae (L675, L944, LI 103, 
LI 111 and LI 246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large 
Array. These objects were previously observed with the AMI Small Array to have an excess 
of emission at microwave frequencies relative to lower frequency radio data. In L675 we find 
a flat spectrum compact radio counterpart to the 850 /J.m emission seen with SCUBA and 
suggest that it is cm-wave emission from a previously unknown deeply embedded young pro- 
tostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 fim emission 
seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that 
it arises from a population of PAH molecules, which also give rise to the cm-wave emission 
through spinning dust emission. 

Key words: Radiation mechanisms: general - ISM: general - ISM: clouds - stars: formation 



1 INTRODUCTION 

The complete characterization of microwave emission from spin- 
ning dust grains is a key question in both astrophysics and cosmol- 
ogy. It probes a region of the electromagnetic spectrum where a 
number of different astrophysical disciplines overlap. It is impor- 
tant for CMB observations in order to correctly characterise the 
contaminating foreground emission; for star and planetary forma- 
tion it is important because it potentially probes a regime of grain 
sizes that is not otherwise easily observable. 

Although a number of objects have now been found to exhibit 
anomalous microwave emission, attributed to spinning dust, it is 
still unclear what differentiates those objects from the many other 
seemingly similar targets that do not show the excess. In the spe- 
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cific case of dark clouds the recent AMI sample (AMI Consortium: 
Scaife et al. 2009; hereinafter Paper I) of fourteen Lynds Dark Neb- 
ulae found an excess in only five. 

It has been suggested that cm-wave emission from spinning 
dust is emitted by a population of ultra-small grains (Draine & 
Lazarian 1998). These ultra-small grains are thought to exist mainly 
in the form of single polycyclic aromatic hydrocarbon (PAH) 
molecules. PAH molecules are generally detected through their nar- 
row line emission features in the MIR. For these emission features 
to be observed the PAH molecules must be exposed to a strong 
source of UV flux. Since this flux is generally absent in the case 
of dark clouds, the microwave emission from the rotation of PAH 
molecules may be the only way to study the very small grain pop- 
ulation in these objects. 

It is also known that radio continuum emission in dark clouds 
may arise from ionized gas associated with a stellar outflow. When 
a luminous star is present this arises either as the result of a com- 
pact Hll region or an ionized stellar wind. In the case of very young 
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low luminosity stars radio continuum emission may be also be de- 
tected. In this instance it is generally attributed to the presence of 
a partially ionized (0.02 < x e 0.35; Bacciotti & Eisloffel 1999) 
stellar wind (Wright & Barlow 1975; Panagia & Felli 1975), or 
possibly a neutral wind which has been shock-ionized further from 
the central source by impacting on a dense obstacle (Curiel et al. 
1989). 

In this paper we present follow-up observations of the five 
AMI Small Array (SA) spinning dust detections (Paper I) at higher 
resolution with the AMI Large Array (LA) over the same frequency 
range. All co-ordinates in this paper are J2000.0. 



Table 1. AMI LA Lynds Dark Nebulae. Column [1] Name of cloud, [2] 
Right Ascension, [3] Declination, [4] AMI LA synthesized beam FWHM 
major axis, [5] AMI LA synthesized beam FWHN minor axis, and [6] r.m.s. 
noise fluctuations on the combined channel map. 



Name 


RA 


Dec 


A0 m aj 


A0 m in 


Grms 




(J2000) 


(12000) 


(arcsec) 


(arcsec) 


bm > 



L675 19 23 52.6 11 07 39 49.9 27.4 35 

L944 21 17 40.8 43 18 08 36.5 31.2 31 

LI 103 21 42 10.2 56 43 44 32.0 26.3 25 

Lllll 21 4027.1 5748 10 39.4 30.8 29 

L1246 23 25 30.1 63 38 30 31.2 26.9 25 



2 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(SA) and one of eight 13 m antennas (LA), both sited at Lord's 
Bridge, Cambridge (AMI Consortium: Zwart et al. 2008). The 
telescope observes in the band 13.5-17.9 GHz from which eight 
0.75 GHz bandwidth channels are synthesized. In practice, the two 
lowest frequency channels (1 & 2) are not generally used due to a 
lower response in this frequency range and interference from geo- 
stationary satellites. 

Observations of five Lynds dark nebulae selected from the 
original AMI SA sample were made in 2009 February-March us- 
ing the AMI LA. The co-ordinates of these fields are listed in Ta- 
ble 1 along with the size of the AMI LA synthesized beam towards 
each object and the r.m.s. noise measured outside the primary beam 
on the CLEANed maps. We note that the AMI LA observation of 
LI 246 is towards the north-east of this cloud where anomalous 
emission was detected by the AMI SA and does not cover the same 
area as the original SCUBA observation. 

Data reduction was performed using the local software tool 
REDUCE, see Paper I for more details. Flux calibration was per- 
formed using short observations of 3C286 near the beginning and 
end of each run. We assumed I+Q flux densities for this source 
in the AMI LA channels consistent with the frequency dependent 
model of Baars et al. (1977), ~ 3.3 ly at 16 GHz. As Baars et al. 
measure I and AMI LA measures I+Q, these flux densities include 
corrections for the polarisation of the calibrator source derived by 
interpolating from VLA measurements. A correction is also made 
for the changing intervening air mass over the observation. The 
phase was calibrated using interleaved observations of calibrators 
selected from the lodrell Bank VLA Survey (JVAS; Patnaik et al. 
1992). After calibration, the phase is generally stable to 5° for chan- 
nels 4-7, and 10° for channels 3 and 8. The FWHM of the primary 
beam of the AMI LA is fa 6' at 16 GHz. 

Reduced data were imaged using the AIPS data package. 
CLEAN deconvolution was performed using the task IMAGR which 
applies a differential primary beam correction to the individual fre- 
quency channels to produce the combined frequency image. De- 
convolved maps were made from both the combined channel set, 
see Fig. 1, and for individual channels. The broad spectral cover- 
age of AMI allows a representation of the spectrum between 14.3 
and 17.9 GHz to be made independently of other telescopes and 
in what follows we use the convention: S v~ a , where S is flux 
density, v is frequency and a is the spectral index. All errors are 
quoted to 1 a. 
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Figure 2. L675 Source A: data points are flux densities from AMI LA chan- 
nels 3-8. The best-fit spectral index of a = 0. 10 ± 0.36 is shown as a dashed 
line. 

3 RESULTS 

L675: The AMI LA observations of L675 show two obvious re- 
gions of compact emission, see Fig. 1. The first of these, slightly 
offset from the pointing centre, is coincident with both the peak 
of the AMI SA emission and also the compact emission seen at 
850 iim by the SCUBA instrument (Visser et al. 2001; 2002). We 
denote this source "A" (19 h 23 m 50?5, +11°07'44"). The second, 
just outside the LA primary beam FWHM to the north-east, is co- 
incident with the probable extragalactic point source identified as 
"B" (19 h 24 m 02?4, +11°10'54") in the original AMI SA observa- 
tions (Paper I). 

Source A is completely unresolved by the AMI LA and shows 
a flat spectrum across the AMI band, aJJ/f = 0.10 ±0.36, consis- 
tent with free-free emission, see Fig. 2. This spectral index differs 
considerably from that measured by AMI SA. This is because the 
LA is not sensitive to the large scale emission seen with the SA. In- 
deed it seems likely that the emission seen by the two arrays arises 
from completely different sources. 

L944: The original AMI SA observations of L944 revealed 
a compact region of emission to the north of the cloud, coincident 
with one side of the protostellar outflow. AMI LA observations, see 
Fig. 1, reveal this emission arises not from a point-like object but 
rather from a diffuse region of emission, the peak of which occurs 
at 21 h 17 m 42?6, +43°20'08". We estimate the flux spectrum by in- 
tegrating the flux density from the primary beam corrected channel 
maps within a two arcminute radius of the LA pointing centre. This 
shows a steeply rising spectrum with ajjf = —2.1 ±0.5. This is 
consistent with that found from the AMI S A data, however this cor- 
respondence is not meaningful as the low signal to noise in the SA 
data precludes a precise estimate. The flux density found towards 
this region in the LA map is only marginally lower than that found 
from the comparatively coarser resolution SA map. This implies 
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L675 L944 




23 h 26 m 00 s 25 m 45 s 25 m 30 s 25 m 15 s 25 m 00 s 



Figure 1. AMI LA combined channel data is shown as grayscale in units of mJy/beam, grey conrtours at —6, —3, ±2,±lCT and black contours at 3, 6, 12, 24 G 
etc. SCUBA 850 jlm data is shown as red contours with levels as in Visser et al. (2001; 2002) for the all clouds except L1246. The AMI LA observation of 
LI 246 does not cover the region observed by SCUBA. AMI SA data is shown as blue contours with levels as in Paper I. The AMI LA primary beam FWHM 
is shown as a circle and the synthesized beam as a filled ellipse in the bottom left corner. 
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Table 2. Integrated flux densities in mJy for AMI LA observations of L675, L944 and L1246. Errors are calculated as 5 = \/ (0.055) 2 + o} m!i , where o~ nTls is 
the r.m.s. noise in the individual channel map. 













Freq. (GHz) 












Name 


14.3 


15.0 




15.7 




16.4 




17.2 




17.9 


a 


L675 
L944 
LI 246 


2.74 


±0.34 2.24 
2.03 
0.62 


±0.12 
±0.11 
±0.18 


2.26 
2.40 
0.55 


±0.12 
±0.13 
±0.03 


2.12 
2.89 
0.62 


±0.11 
±0.15 
±0.03 


2.42 
2.86 
0.58 


±0.14 
±0.15 
±0.04 


2.22 ±0.12 


±0.10±0.36 
-2.11 ±0.49 
-0.40 ±0.82 



that the emission comes not from one smooth extended region that 
is partially resolved out by the LA baselines, but from a collection 
of smaller fragments or filaments. These fragments are unresolved 
by either array, although the granularity becomes more evident in 
the higher frequency channels of the LA. The amount of flux lost 
in channels 5 to 8 relative to channel 4 is significantly smaller than 
would be expected from a Gaussian source of similar dimensions. 

L1103 and Lllll: AMI LA observations of LI 103 and 
LI 1 1 1 do not show any distinct regions of compact cm- wave emis- 
sion. The diffuse patches of low level emission present within the 
primary beam towards both sources are indicative of larger scale 
structures which have been resolved out by the synthesized beam. 
We can provide an estimate of the flux density seen towards these 
objects with the LA by fitting and removing a tilted plane baselevel 
at the primary beam FWHM. From the combined channel data this 
gives Si6 = 5.1 ±0.6mJy and S\(, = 2.4±0.3mJy for LI 103 and 
Lllll, respectively. These values indicate that the flux loss is con- 
siderable: approximately 45% and 96%. The sensitivity of these 
LA observations is much greater than those of the SA and it is 
possible that the patchy emission in these fields corresponds to en- 
hancements in the extended emission which are below the detection 
threshold in the SA data. 

L1246: AMI SA observations towards L1246 did not show 
any excess emission coincident with the SCUBA identification of 
the dark cloud, but did reveal a region of emission ~ 2' to the north- 
east of the cloud, in a region not covered by the SCUBA map, 
which had no counterpart in the lower frequency observations. AMI 
LA observations of this NE region show an arc of emission (peak: 
23 h 25 m 20?4, ±63°38'40"), see Fig. 1. We assess the spectral be- 
haviour of this object in two ways. Firstly, we estimate the flux 
density of the arc itself. We fit and remove a tilted plane baselevel 
within an irregular polygon drawn around the object and integrate 
the remaining flux. Secondly, we fit a tilted plane baselevel to a 
circle at the primary beam FWHM and integrate all the flux above 
this baselevel within that radius. Both methods give consistent re- 
sults, as might be expected since the primary beam is relatively 
empty otherwise. From the first method we find a spectral index 
a\H = -0.40 ±0.87, and from the second a>^ = -0.47 ±0.82. 



4 DISCUSSION AND CONCLUSIONS 

L675 and LI 246 have archival Spitzer IRAC data, which shows in 
both cases a significant amount of emission in Band 4 (6.4-9.4 flm) 
and very little in the other three (3.2-3.9, 4.0-5.0 and 4.9-6.4 flm, 
respectively). In the case of L675 this emission is present on a very 
large scale, see Fig. 3. The emission seen at 16 GHz with the AMI 

5 A appears on a similar scale, however the small field of view of the 
Spitzer data precludes a more detailed comparison. LI 246 shows 
an arc of emission at 16 GHz which is also evident in Spitzer IRAC 
Band 4, see Fig. 4. This emission is again not present in Bands 1-3. 




19 h 24 m 30 s 24 m 15 s 24 m 00 s 23 m 45 s 23 m 30 s 23 m 15 s 

Figure 3. L675: AMI LA combined channel data is shown as white contours 
at 3, 6, 12 o~ etc. Spitzer Band 4 data is shown as greyscale in MJy/sr, and 
is saturated at both ends to emphasise the structure present. AMI SA data is 
shown as blue contours as in Fig. 1 . 




23 h 26 m 00^5 m 45 s 25 m 30 s 25 m 15 s 25 m 00 s 



Figure 4. L1246: AMI LA combined channel data is shown as white con- 
tours at 1, 2, 4, 8 a etc. Spitzer Band 4 data is shown as greyscale in MJy/sr, 
saturated at both ends of the scale to emphasise the structure present. The 
AMI LA primary beam is shown as a circle and the synthesized beam as a 
filled ellipse in the bottom left corner. 

In Band 4 it is present as an arc, coincident with that seen at 16 GHz 
in the AMI LA data. 

Spitzer Band 4 contains two of the PAH emission lines, in- 
cluding the strongest (7.7 flm). Of the three other Spitzer bands 
only Band 1 contains an emission line (3.3 flm) and for ionized 
PAHs this line is expected to be significantly weaker. It is probable 
therefore that the MIR correlated cm-wave data seen in the AMI 
maps is a consequence of spinning dust emission from a popula- 
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tion of ionized PAH molecules. Neutral PAH molecules do not in 
general possess a permanent dipole moment and are therefore not 
expected to have rotational emission (Tielens 2008). This emission, 
the mechanism of which is described in detail by Draine & Lazar- 
ian (1998), arises from the intrinsic dipole moments of small dust 
grains, most likely to be PAH molecules, which emit power when 
they rotate. This rotation has a variety of contributing factors, the 
relative importance of which varies with grain environment. How- 
ever, in the majority of cases excitation through collision with ions 
is predominant. 

In the case of L675 A, we must consider the possibility that we 
are observing a coincidental extragalactic radio source. Using the 
extended 9C survey 15 GHz source counts (Waldram et al. 2009), 
wheren(S) = 51(S/Jy)~ 2 ' 15 Jy~' sr~ 1 , the probability that a source 
with flux density greater than 2 mjy lies within the FWHM of the 
AMI LA primary beam is 0.12, and only 0.01 within the SCUBA 
field. It is likely therefore that the radio source L675A is associated 
with the SCUBA core. 

A further question is whether the cm-wave emission might be 
explained by thermal (Planckian) dust emission. A single greybody 
spectrum with a dust temperature, Ri 27 K, might be used to ex- 
plain the LA flux density, however it would require a /3 of 0.6. Such 
a value would be unusual even for objects known to possess flat- 
tened dust tails, such as protoplanetary disks. This simple fit also 
neglects the flux lost by the AMI LA baseline distribution. SA ob- 
servations have already shown this source to possess a significant 
amount of extended emission which would make this scenario even 
more unlikely. 

The presence of a neutral or partially ionized wind from an 
outflow source that has been shocked through encountering a dense 
obstacle (Torrelles et al. 1985; Rodriguez et al. 1986) is used to 
understand the spectral indices seen towards exciting sources in 
the radio regime (Curiel et al. 1990; Cabrit & Bertout 1992). This 
model allows a spectral index range of 0.1 (optically thin) to -2 
(optically thick), which explains results which deviate from the 
value of a = —0.6 required by a spherically symmetric ionized 
wind (Wright & Barlow 1975; Panagia & Felli 1975). Using this 
model as described in Curiel et al. (1989; 1990) the radio emis- 
sion is expected to be optically thin (t = 0.1), consistent with the 
spectral index seen across the AMI band. Assuming a distance of 
300 pc and a stellar wind with a wind speed of 200 km s _1 , we can 
calculate that the AMI flux densities towards L675A are consistent 
with a mass loss of 3.5 x 10~ 7 M Q yr~' . A mass loss such as this 
implies a mechanical luminosity from the wind of L mec h fa 1.1 L©, 
comparable to the values found by Curiel et al for L1448. 

The nature of the emission seen towards L944 with the AMI 
LA is uncertain. The spectral index of this emission is consis- 
tent with spinning dust emission or alternatively the optically thick 
component of free-free spectrum. Such a free-free spectrum might 
be exhibited at 16 GHz by ultra-compact Hll regions. However a 
turn-over frequency above 16 GHz would have an extremely high 
mass and should therefore be obvious in sub-mm observations. This 
needs to be confirmed by either higher radio frequency measure- 
ments in order to measure the optically thin region of the spectrum 
and the turn-over, or sub-mm measurements to place constraints on 
the mass of such a region. 

In conclusion, we have used the AMI LA to observe a sample 
of five Lynds Dark Nebulae selected as candidates for spinning dust 
emission from the AMI SA sample of Lynds Dark Nebulae (Paper 
I). Towards two of these clouds (LI 103 and LI 1 1 1) we detect only 
patchy diffuse emission characteristic of the presence of a larger 
structure which has been mostly resolved out. 



Towards L675 we have observed flat spectrum compact cm- 
wave emission coincident with the SCUBA 850 flm emission from 
the same region. These characteristics suggest that this source is 
associated with a stellar wind from a deeply embedded young pro- 
tostar. 

We detect extended cm-wave emission to the North of the 
L944 SMM-1 protostar which displays spectral behaviour consis- 
tent with either spinning dust, or alternatively a collection of ultra- 
compact Hll regions. 

L1246 shows an arc of cm-wave emission which is coincident 
with emission seen in Spitzer Band 4. We suggest that this is an 
example of emission from a population of PAH molecules, seen in 
emission lines in the Spitzer data, and emission as a consequence 
of rapid rotation of the molecules in the cm-wave data. 
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